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RESUMEN 

Se presentan dos modelos simples de columna para representar la incidencia del efecto de 
corte y de flexion y la interaccion entre estos y la carga normal en la respuesta global de columnas 
solicitadas por estos esfuerzos de corte, flexion y axial. Los modelos constan de barras elasticas con 
resortes de flexion no lineales en los extremos y resortes no lineales de corte, concentrados en el 
centro de la luz para uno de ellos y en los extremos para el segundo modelo. Todos los resortes 
siguen modelos histereticos que representan los efectos principales de flexion y corte y con 
envolventes obtenidas en forma simple por metodos usuales para estructuras de hormigon armado El 
acoplamiento entre la flexion y corte se consigue trabajando con la matriz de rigidez y la interaccion 
con la fuerza axial se obtiene modificando las envolventes resultante de la variacion de la carga axial. 
Se puede concluir que los modelos simples representan muy bien la interaccion entre los esfuerzos 
de flexion y corte y la degradacion de resistencia con los ciclos de carga/descarga y resulta de facil 
implementacion en programas de analisis no lineal. 



ABSTRACT 

Two simple models of column are presented, to represent the effects of shear and bending 
forces and the interaction between them and the axial forces. The models are elastic bars with 
nonlinear flexural springs at both ends and nonlinear shear springs at the ends for one of the models 
and one central nonlinear shear spring in the other case. All the springs follow hysteretic rules which 
represent the principal effects of flexure and shear, with envelopes which are obtained with simple 
procedures for reinforced concrete structures. The coupling of bending and shear is given in the 
stiffness matrix and the interaction with the axial force is given by the modification of the envelopes 
due to the variation of it. The models represent, in an accurate way, the coupling between bending and 
shear forces and the interaction of them with the axial force, also the strength degradation due to the 
unload/reload cycles is well represented. The models are simple to implement in a nonlinear analysis 
program. 
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1. INTRODUCCION 

Las fallas de corte en elementos de hormigon armado son del tipo fragil. Si bien actualmente 
se trata de evitar este tipo de fallas usando procedimientos de diseno por capacidad, existen 
estructuras que han sido disenadas con anterioridad a la aplicacion de estos criterios o bien, debido a 
que se consideraba un comportamiento diferente, se le daba poca importancia estos esfuerzos. Esto 
ultimo ocurre en las estructuras de mamposteria encadenada sujetas a acciones horizontales, donde 
si bien los encadenados, especialmente los verticales, estan sometidos principalmente a esfuerzos 
axiales, los esfuerzos de flexion y corte tienen, en algunos casos, valores elevados (Combescure 
2006, Torrisi y Crisafulli 2010). 

Existen en la actualidad diversos modelos para representar el comportamiento de columnas 
de hormigon armado, que van desde modelos muy simples hasta algunos muy complejos en donde 
la cantidad de datos a considerar es muy elevado (Ceresa et al. 2009). Algunos de estos modelos 
son los llamados modelos de fibras, los cuales subdividen a la seccion en una cierta cantidad de 
fibras de hormigon y acero cada una de ellas con su correspondiente modelo histeretico, ademas, 
subdividen al elemento en varias secciones para su posterior integracion, donde cada una de ellas es 
representada por la seccion de fibras. Un programa muy difundido que utiliza este metodo es 
SEISMOSTRUCT (Seismosoft 2007). Si bien estos modelos representan en forma adecuada el 
comportamiento flexional del elemento, algunos son incapaces de representar el efecto del corte . 
Por tal motivo se han desarrollado nuevos modelos (Ceresa et al. 2009) que consideran estos 
efectos, sin embargo, el costo computacional y la carga de datos hace, que su uso, sean muy 
complejos.. 

Algunos autores (Filippou y Issa 1988, Scapone et al. 1992 han propuesto elementos mas 
simples en donde realizan una mezcla entre modelos de fibras y el comportamiento global del 
elemento, definido por un ciclo histeretico y una envolvente de esfuerzos, sin embargo, persiste el 
problema de la definicion de la seccion del elemento para utilizar el modelo de fibras. 

En este trabajo se presentan dos modelos simple que tiene en cuenta los esfuerzos de 
flexion y de corte y una interaccion entre ellas y las fuerzas axiales. La formulacion de los modelos 
esta dada en terminos de flexibilidad y el comportamiento flexional y de corte se realiza mediante la 
definicion de las envolventes momento-curvatura y corte-distorsion que siguen leyes histereticas 
previamente definidas. 



2. FORMULACION DE LOS ELEMENTOS 



Para modelar el comportamiento de los elementos de hormigon armado sometidos a 
esfuerzos de corte y flexion acoplados se han desarrollado 2 modelos simples. El primero de ellos, 
elemento EMV1a, es una viga de Euler-Bernoulli con energia de corte donde el comportamiento de 
flexion y corte esta dado por resortes no lineales en los extremos. El segundo elemento, EMV1b, es 
tambien una viga de Euler-Bernoulli con energia de corte en donde la flexion se modela mediante 
resortes no lineales en los extremos y el cote mediante un resorte central, tambien con 
comportamiento no lineal. 

Se plantean dos modelos para evaluar las ventajas y desventajas en cada uno de ellos. En el 
primero de ellos, la ventaja principal es la simple implementacion. En este caso, la matriz de 
flexibilidad total del elemento es la suma directa de matrices de flexibilidad para los distintos efectos., 
sin embargo, al considerar el esfuerzo de corte distribuido en el elemento (representado por dos 
resortes extremos), la degradacion se produce en todo el elemento y no se pueden conocer ciertas 
variables como por ejemplo los desplazamientos de corte , ademas no representa en forma fisica la 
fisura de corte. El segundo modelo tiene la ventaja de poseer un resorte fisico que permite conocer 
el desplazamiento de las secciones y representa la fisura de corte en forma mas exacta. Ademas, 
este modelo degrada su resistencia y rigidez a corte solo en la zona del resorte. Sin embargo, la 
implementacion es mas compleja ya que requiere una condensacion de los grados de libertad 
internos del elemento (correspondientes al elemento resorte de corte), a los grados de libertad 
externos del elemento. 
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La figura 1 muestra esquematicamente los elementos EMV1a y EMV1b. En primer lugar se 
muestra el elemento EMV1a con sus grados de libertad y las flexibilidades de los resortes de flexion 
fsi y f S 2, y de corte f c i. En segundo lugar se presenta el elemento EMV1b previo a ser condensados 
los grados de libertad internos y finalmente el elemento EMV1 b una vez realizada la condensacion de 
los grados de libertad internos. En esta figura tambien se muestran los resortes de flexion y corte con 
sus flexiblidades f s i y f S 2 y la rigidez K v , respectivamente. 
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Figura 1 Elementos EMV1 a y EMV1 b con sus grados de libertad. a) Elemento EMV1 a, b) Elemento 
EMV1 b con los grados de libertad antes de ser condensados. c) Elemento EMV1 b con los grados de 

libertad condensados 



En los elementos mostrados en la figura 1 , los grados de libertad del elemento (1 a 6), son los 
expresados en la ecuacion 1, los cuales estan relacionados con los esfuerzos mostrados en la 
ecuacion 2. 

{u}=[x, y, 9, x 2 y 2 0 2 ] T (1) 
{Q}=[N, V, M, N 2 V 2 Mj (2) 



En la ecuacion 1, x, y, 6, representan los desplazamientos en direccion de la barra (axial, 1 y 
4), perpendicular a la misma (2 y 5) y el giro de los nodos extremos (Flexion, grados de libertad 3 y 6). 



2.1 Elemento EMV1 a 



El primer modelo planteado corresponde a una barra de Euler-Bernoulli con energia de corte. 
La formulacion de este elemento se realiza en terminos de flexibilidad (Flores-Lopez 1993, D'ambrisi 
y Filippou 1999, Leu y Cheng 2000, Al-Haddad 1990, Thonson y Flores-Lopez 2004, Filippou y Issa 
1 998) debido a que resulta mas sencilla la consideracion de los distintos efectos. 

La matriz de flexibilidad total elemental del modelo esta compuesta por 4 submatrices. La 
primera de ellas corresponde a la formulacion pura de Euler-Bernoulli, Feb, obtenida aplicando 
momentos unitarios en los extremos de una viga simplemente apoyada. La segunda tiene en cuenta 
la energia de corte, F ec , la cual se obtiene por el corte inducido al aplicar momentos unitarios en los 
extremos de la viga anteriormente mencionada. La tercer submatriz, denominada Fu, incluye la 
contribucion de las rotulas de flexion de los extremos, es decir, la flexibilidad del resorte de flexion y la 
cuarta submatriz, llamada F c ;, es la encargada de aportar el comportamiento no lineal por corte. La 
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ecuacion 3 muestra las distintas matrices de flexibilidad sumadas para obtener la matriz de flexibilidad 
elemental Fdel subelemento EMV1a. 



[Fh[F eb ]+[F ec h[F b i,h[F ci ] 



(3) 



La formulacion de las matrices de flexibilidad elasticas esta dada, por las ecuaciones 4 a 7 
(Aguiar 2003) 
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Donde A es el area de la seccion transversal del elemento, A c el area de corte del mismo, J 
es el valor de momento de inercia de la seccion transversal, Ey G son los modulos de elasticidad y 
de corte del material y U es la longitud del tramo. Los factores f s i y fs2 son las flexibilidades de los 
resortes de flexion y fe la flexibilidad de corte y son definidos segun las ecuaciones 8 y 9 (Carr 2007). 
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En las ecuaciones anteriores n y p, son los factores bilineales para flexion y corte 
respectivamente, definidos como la rigidez tangente instantanea sobre la rigidez inicial. El factor a c 
fue calibrado con ensayos experimentales y tiene en cuenta la variabilidad del modulo de corte y del 
area de corte. Cabe aclarar que el modelo histeretico a utilizar no puede tener comportamiento 
plastico perfecto (Carr 2007), es decir una pendiente igual a cero ya que el factor bilineal seria infinito, 
es por ello que se debe introducir un valor no nulo, positivo o negative 

Una vez obtenida la matriz F se procede a invertir la misma para obtener la matriz K 0 . La 
matriz de rigidez del elemento K e , se obtiene con la ecuacion siguiente: 
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KHsfKIe] no) 

Donde B es la matriz de transformacion definida como: 
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La matriz Ke representa la matriz de rigidez del elemento EMV1a, en coordenadas locales, 
con un orden de desplazmientos segun la ecuacion 1 . 



2.2 Elemento EMV1 b 



El segundo modelo planteado corresponde, al igual que el primero, a una barra de Euler- 
Bernoulli con energia de corte. En este caso se analiza la matriz de flexibilidad de los tramos elasticos 
anterior y posterior al resorte de corte central y la contribucion de los resortes flexionales extremos. 
Estas matrices se invierten para obtener las matrices de rigidez de cada tramo y se ensamblan con la 
matriz de rigidez del resorte colocado a mitad de longitud. La matriz de rigidez resultante posee 12 
grados de libertad, los cuales hay que condensar para llevar a los grados de libertad de los nudos 
extremos (Ascheri et a/ 2009), quedando asi una matriz de rigidez de 6x6. 

La ecuacion 12 expresa la matriz de flexibilidad elastica de cada tramo, de longitud Li, del 
subelemento EMV1b, considerando que el resorte de corte se encuentra a mitad de longitud. 
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Las matrices de flexibilidad que tienen en cuenta la rigidez de los resortes de flexion extremos 
estan dadas en las ecuaciones 13 y 14. 
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La matriz de rigidez elemental de cada tramo Ki y K2, se obtiene invirtiendo las matrices F e i y 
Fe2, definidas como: 



[FeMFehM [F e2 ] = [F e ]+[F b2 ] 

J 



(15) 
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J 

La matriz de rigidez en coordenadas locales de cada tramo se obtiene de la siguiente forma 

K] = [bJ[kM] 

[K e2 ]=[Bj[K 2 )[B,] 
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La matriz de rigidez del resorte central queda expresada como sigue 
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(20) 



donde K y es la rigidez del resorte de corte y K xe y K qe son las rigideces elasticas axial y de flexion 
definidas como sigue 



K xe = 



K 8e = 



2EA 
12EJ 



(21) 



(22) 



La rigidez K y se obtiene como la inversa de la ecuacion 7 para valores de p/distintos de 1. Si 
el valor de p, es igual a 1 , comportamiento elastico, el valor de la rigidez se evalua como V C r/ycr, siendo 
V cr e\ corte de agrietamiento y y C r\a distorsion de corte correspondiente. 

Las matrices K e i, K e 2y K v se ensamblan en serie en una matriz K 0 , colocando la submatriz del 
subelemento K e i, a continuacion la matriz del resorte K v y finalmente la matriz de rigidez 
correspondiente al subelemento K e 2. Una vez realizado el ensamblaje de las submatrices 
mencionadas en la matriz K 0 , se procede a ordenar los grados de libertad de modo que en primer 
lugar se encuentren los grados de libertad correspondientes a los nudos extremos y posteriormente 
los correspondientes a los nodos interiores. Esto se realiza con la ecuacion 23, valiendose de la 
matriz C, definida en la ecuacion 24. 



[K d ] = lc] T [Kjc] 



(23) 



La matriz Kd es una matriz que se encuentra subdivida en 4 submatrices. La primera de 
ellas Km incluye los 6 grados de libertad extremos, Kbb, los 6 grados de libertad internos y Kab y 
Kba que incluyen el acoplamiento entre ambos. 
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[C] 



100000000000 
010000000000 
001000000000 
000000100000 
000000010000 
000000001000 
000000000100 
000000000010 
000000000001 
000100000000 
000010000000 
000001000000 



K] 



[K AA ] [K AB ] 
[K BA ] [K BB ] 



[K e ]=[K M ]-[K AB lK BB V[K BA ] 



(24) 



(25) 
(26) 



Mediante la condensacion, realizada con la expresion 26, se eliminan los grados de libertad 
internos del elemento, aquellos que vinculan el primer tramo de viga con el resorte central y este con 
el segundo tramo de viga. La matriz obtenida por aplicacion de la ecuacion 26 posee los grados de 
libertad ordenados de acuerdo a la expresion 1 . 

Las matrices de rigidez de los elementos K e , definidas por las ecuaciones 10 y 26, se 
ensamblan luego en la matriz de rigidez global de la estructura, previo paso de coordenadas locales a 
coordenadas globales. 



2.3 Solucion incremental y obtencion de deformaciones y fuerzas internas. 



Una vez ensamblada las matrices de rigidez de los elementos en el sistema, se procede a 
resolver el sistema de ecuaciones. Esta solucion se obtiene en forma incremental, obteniendo 
incremento de desplazamientos para los nodos del macroelementos, Au m , y por consiguiente 
incrementos de desplazamientos para cada subelemento, Au e - Ademas, si la solucion es del tipo 
implicita, -se puede utilizar algun procedimiento tipo Newton-Raphson (Filippou et al 1992), para 
verificar el equilibrio global y local. Alternativamente se puede utilizar un procedimiento explicito para 
resolver el sistema y los elementos. 

El incremento de fuerzas internas, AQ, viene dado por la ecuacion 27, donde Au e esta 
definido en la ecuacion 28. 

{AQ}=[AA/ 1 A^ AM, AN 2 AV 2 AMj =[K e ]{Au e } (27) 
{Au e }=[Ax, Ay 1 AG, Ax 2 Ay 2 A0 2 J (28) 



A partir del incremento de fuerzas internas AQ, se calcula el incremento de deformaciones 
internas (Carr 2007), que para el caso de curvaturas se obtiene como se indica a continuacion. 
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* j. A/W 1 

A ^ = ^7 <29) 

A0 2 = — 2 (30) 
r 2 EJ 

Para el caso particular del elemento EMV1a, el incremento de la distorsion de corte se calcula 
de acuerdo con la ecuacion 31 . 

AV< 

= ^ (31) 

Si se trabaja con el elemento EMV1b, el incremento de la distorsion de corte se obtiene como 
se explica a continuacion. Primero se obtiene la matriz de desplazamientos internos K, (ecuacion 32), 
luego se calcula el desplazamiento en el resorte central, segun la ecuacion 33 (Ascheri et al 2009) y 
posteriormente se evalua el esfuerzo en dicho resorte con la ecuacion 34, para finalmente obtener el 
valor buscado del incremento de distorsion de corte, ecuacion 31 . 

[K,] = -[KbbY[K ba ] (32) 

M=[K,^U g } (33) 

{AQ}=[KjAu r } (34) 



Con los valores incrementales de las deformaciones internas se entra en los diagramas 
histereticos correspondientes y se obtienen los valores de rigidez tangente instantanea, factores 
bilineales y fuerzas internas. Si dichas fuerzas no son iguales, dentro de un margen de error 
establecido, a las externas, se procede a realizar iteraciones de equilibrio en el incremento para 
asegurar la convergencia. 

Para tener en cuenta la interaccion entre las fuerzas axiales y de corte y las fuerzas axiales y 
de flexion, se modifican las envolventes de dichos esfuerzos, (corte y flexion), con el valor de la 
fuerza axial en el elemento. Es por ello que se ingresa primero al modelo histereticos para la fuerza 
normal y luego con ese valor se modifican las envolventes y se ingresa con los incrementos de 
curvatura y distorsion de corte. 



3. COMPORTAMIENTO DEL MODELO Y LEYES HISTERETICAS 



El comportamiento de cada resorte no lineal esta caracterizado por su ciclo histeretico, el cual 
relaciona la deformacion interna (curvatura, distorsion de corte o deformacion axial) con el esfuerzo 
caracteristico (momento, corte o fuerza axial). Estos ciclos histereticos son definidos en funcion del 
comportamiento observado para cada esfuerzo. En los elementos propuestos se utiliza el modelo de 
Takeda (figura 2) para representar el comportamiento de los resortes de flexion, y un modelo trilineal 
con estrangulamiento de lazo (pinching), mostrado en la figura 3, para representar el comportamiento 
bajo esfuerzos de corte. El comportamiento axial se modela segun una ley histeretica bilineal (figura 
4) en traccion, con rigidez postfluencia r/K 0 y un comportamiento lineal en compresion. Como 
parametro comun, en la definicion de todos los modelos histereticos se ha considerado una 
degradacion de rigidez en la descarga, modificando la rigidez inicial por un factor menor que la 
unidad. 
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Figura 4 Modelo bilineal para los esfuerzos axiales 
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A partir de los modelos histereticos correspondientes a cada uno de los resortes, es posible el 
conocer la rigidez instantanea para un estado de deformacion dado. La rigidez en cada instante se 
define mediante los factores n, para flexion, p ( para corte y ri para la carga axial, segun muestra la 
figura 5. Con estos factores se puede modifican las matrices de flexibilidad y rigidez de los resortes. 

Con la rigidez instantanea y la definicion de la flexibilidad de los resortes es posible conocer 
la matriz de flexibilidad y por lo tanto la matriz de rigidez, de cada elemento en cada paso de carga, 
como se explico previamente . 

Para considerar la variacion de las propiedades de flexion con la carga axial se construye el 
diagrama de interaccion de flexion-axial, definido al menos por 4 puntos (Park y Paulay 1975, Llopiz 
2007). Para las propiedades de corte se puede realizar algo similar, modificando las resistencias en 
funcion de la carga axial existente, utilizando por ejemplo, la formulacion del USCD (Calvi et a/ 2005). 
La figuras 6a y 6b muestran la respuesta en flexion y corte cuando se considera la carga axial 
variable y cuando la misma se mantiene constante. En ambos casos se ha considerado que cuando 
las deformaciones son positivas la carga axial aumenta su valor de compresion y en forma similar, 
cuando las deformaciones son negativas, se considera que la carga axial aumenta su valor de 
traccion. 




-60,0 - 

-80,0 - 
Curvatura (rad/m) 
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-60,0 - 



-80,0 J 
Distorsion (mm/mm) 



Figura 6 Respuesta histeretica de flexion a) y de corte (b, considerando la carga axial constante y 

variable 



Ademas de definir el diagrama de interaccion para flexion y corte es necesario modificar los 
ciclos histereticos para evitar un salto brusco en rigidez ante una variacion de la carga axial gradual. 
Se puede apreciar que el modelo capta la variacion de la carga axial sin cambios bruscos (saltos de 
carga) en el comportamiento, sino que se modifica la rigidez instantanea con dicha variacion de 
carga. 



4. VALIDACION DEL MODELO 



El modelo de columna desarrollado se ha validado en forma monotonica y ciclica con 
resultados experimentales de pilas de puentes (Calvi et a/ 2005) y con paneles sandwich M2 (Bounas 
y Pavese 2010), todos ensayados en los laboratorios del EUCENTRE, en Pavia. Ademas se comparo 
con los resultados experimentales de columnas ensayadas por Sezen, en el aho 2002. Los distintos 
modelos han sido sometidos a una historia de desplazamientos cfclicos crecientes e independientes 
del tiempo, impuestos en la cima. 

Las caracteristicas geometricas de las columnas de puente, ensayadas por Calvi et al. (2005) y 
Peloso (2003), se pueden ver en las figuras 7 y 8 y la tabla 1 muestra un resumen de las 
propiedades. 
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Figura 7 Secciones transversales de las columnas S (tipo 1 ) y T (tipo 2). 
(Calvi etal. 2005, Peloso 2003) 




Propiedades geometricas de los especiimenes S y T. 
Figura 8 Columnas S y T. (Calvi et al. 2005) 
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Tabla 1 Propiedades de las columnas de puente de Calvi et al (2005). 



Columna 


f'c (MPa) 


Altura 

(mm) 


Refuerzo de 
la seccion 


Carga axial 
(kN) 


Modo de 
colapso 


S250 


35.0 


900 


Tipo 1 
pu=1.07% 
p s i=0.13% 


250 (v=0.06) 


Corte 


S500 


23.7 


500 (v=0.19) 


T250 


30.3 


1350 


Tipo 2 
pu=1.77% 
p s i=0.25% 


250 (v=0.07) 


Corte 
combinado 
con flexion 


T500 


29.7 


500 (v=0.15) 



En las figuras 9 se observa la respuesta ciclica obtenida analiticamente con los elementos 
EMV1a y EMV1b comparada con la respuesta experimental de las pilas de puentes ensayadas por 
Calvi et al (2005). Se puede ver que en todos los casos los modelos representan en forma adecuada 
la degradacion de rigidez y resistencia vista experimentalmente. Ademas se nota el efecto de 
pinching, el cual es tfpico en elementos dominados por esfuerzos de corte. 



La figura 10 muestra un corte de los paneles M2 ensayados en EUCENTRE (Bournas y 
Pavese 2010). Estos paneles estaban compuestos por un centro de espuma de poliestireno sobre la 
cual se colocaba una malla de acero en ambos costados yr posteriormente se recubrfa con hormigon. 
Este sistema mostro un comportamiento donde los lazos sufrian un gran estrechamiento en los ciclos 
(pinching), tfpico de elementos dominados por esfuerzos de corte. La envolvente de corte utilizada en 
estos analisis fue obtenida en forma experimental, a partir de los valores de tensiones de corte y 
modulos de elasticidad transversal calibrados en trabajos previos. (Bournas y Pavese 2010) 
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Figura 9 Comparacion con resultados experimentales - respuesta ciclica. 
a)S250-EMV1a, b)S250-EMV1b, c)S500-EMV1a, d)S500-EMV1b, e)T250-EMV1a, f)T250-EMV1b, 

g)T500-EMV1a, h)T500-EMV1b 
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Figura 10 Descripcion panel M2 (Bournas y Pavese 2010) 
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Tabla 2 Propiedades de los paneles M2 (Bournas y Pavese 201 0) 



Especimen 


Dimensiones (m) 


Dimensiones de 
abertura (m) 


Carga axial 
(kN) 


P3-150 


3.00x3.00 




150 


P3-300 


3.00x3.00 




300 


P4-150 


4.00x3.00 




150 


P4-300 


4.00x3.00 




300 


P3-150-W 


3.00x3.00 


Ventana (1.2x1.20) 


150 


P3-300-W 


3.00x3.00 


Ventana (1.2x1.20) 


300 


P3-150-D 


3.00x3.00 


Puerta (2.1 x0.9) 


150 


P3-300-D 


3.00x3.00 


Puerta (2. 1 x 0.9) 


300 



Las figuras 11 muestran las envolventes monotonicas obtenidas para los paneles M2 con el 
elemento EMV1a. Se observa en todos los casos que la respuesta analitica es una envolvente de la 
respuesta experimental, y es capaz de tomar la degradacion de resistencia que se genera en los 
paneles. 
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e) f) 




Desplazamiento (mm) Desplazamiento (mm) 



Figura 11 Comparacion con resultados experimentales - monotonico. 
a)P3-150, b)P3-300, c)P3-150-D,d)P3-300-D, e)P3-150-W, f)P3-300-W, g)P4-150, h)P4-300 

En 2002, Sezen ensayo una serie de columnas de hormigon armado con el objeto de estudiar 
el comportamiento sismico de las mismas. El esquema general de las columnas ensayadas se puede 
ver en la figura 12, donde se muestra una seccion transversal tipica y una elevacion de la columna. 
Algunas de estas columnas se modelaron con los elementos propuestos y se calculo la respuesta 
monotonica de las mismas. Esta respuesta se presenta en la figura 13 y muestra una muy buena 
correlacion con las envolventes experimentales obtenidas. 
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Igual refuerzo al de la viga inferior 



5" 9 barras longitudinales 



W<0 9 B" 




tubos de 2" de 
tliametro 



estribos #3@4" O.C. 
□ 

estribos #3@6" O.C. 



tubos de 2" de 
darnetro 



Figura 12 Columnas ensayadas por Sezen (2002) 

En las figuras 13 se observa la respuesta estatica envolvente de las columnas ensayadas por 
Sezen (2002). Se ve que el modelo propuesto es capaz de representar adecuadamente el 
comportamiento no lineal, por corte y flexion de estos elementos. De esta manera, el elemento 
propuesto es apto para modelar estructuras de hormigon armado en donde se esperen fallas por 
flexion, corte, o combinacion de ambos esfuerzos. 
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b) 



:Experimental 




10 I 1 1 1 1 1 1 

-60 -40 -20 0 20 40 60 

Desplazamiento (mm) 



c) 



300 




10 I 1 1 1 1 1 1 1 

-60 -40 -20 0 20 40 60 80 



Desplazamiento (mm) 
Figura 13 Comparacion de resultados de las columnas ensayadas por Sezen (2002) 



La comparacion con resultados experimentales indica que los modelos propuestos son aptos 
para modelar estructuras de hormigon armado donde se esperan fallas por flexion, corte, o 
combinacion de ambas. Ademas el elemento permite la interaccion de los esfuerzos de corte y 
flexion con las fuerzas axiales, situacion que los hace adecuados para representar el efecto de 
cargas axiales variables y la variacion de la resistencia con las mismas. 
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5. CONCLUSIONES 



Se presentan en este trabajo dos elementos columna propuestos para analizar estructuras de 
hormigon armado sometidos a fuerzas axiales, de corte y de flexion. La principal ventaja de estos 
modelos es que permiten considerar la interaccion entre los esfuerzos axiales y de corte y entre los 
esfuerzos normales y de flexion. El primero de ellos, denominado EMV1a, incluye resortes no 
lineales de flexion y de corte en los extremos, y el segundo con dos resortes de flexion en los 
extremos y un resorte central de corte, llamado EMV1b. La diferencia fundamental entre ambos es la 
simplicidad de implementacion del primer modelo y la capacidad de representar fisicamente el 
desplazamiento de la fisura con una mayor dificultad de implementacion del segundo modelo. 

El comportamiento de cada resorte se representa mediante leyes histereticos que siguen una 
envolvente obtenida con los metodos comunes de analisis de secciones de hormigon armado, las 
cuales se modifican por efecto de la carga axial. Se ha probado la bondad de los elementos mediante 
comparaciones con ensayos experimentales. Los elementos muestran una adecuada envolvente 
estatica de las respuestas experimentales y es capaz de representar, en forma ciclica, la degradacion 
de rigidez y resistencia vista experimentalmente, con lo cual se concluye que el modelo es adecuado 
para representar el acoplamiento entre las fuerzas de corte y de flexion, con carga axial variable, en 
elementos de hormigon armado. 
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